Static and dynamic evolution of CO adsorption on y-U (1 0 0) surface 


with different levels of Mo doping using DFT and AIMD calculations 


Li Jun-Wei!)*) Jia Wei-min?) Liu Chong!) Lv Sha-Sha?) Wang Jin-tao?) Li Zheng-cao!)* 
1) (Key Lab of Advanced Materials (MOE), School of Materials Science and Engineering, 
Tsinghua University, Beijing 100084, China) 
2) (Xi'an Research Institute of High-Technology, Xi'an 710025, China) 
3) (Key Laboratory of Beam Technology (MOE), College of Nuclear Science and Technology, 
Beijing Normal University, Beijing 100875, China) 


Abstract: Alloys of uranium and molybdenum are considered as the future of nuclear fuel and 
defense materials. However, surface corrosion 18 a fundamental problem in practical applications 
and storage. In this study, the static and dynamic evolution of carbon monoxide (CO) adsorption 
and dissociation on y-U (1 0 0) surface with different Mo doping levels was investigated based on 
density functional theory and ab initio molecular dynamics. During the static calculation phase, 
parameters, such as adsorption energy, configuration, and Bader charge, were evaluated at all 
adsorption sites. Furthermore, the time-dependent behavior of CO molecule adsorption were 
investigated at the most favorable sites. The minimum energy paths for CO molecular dissociation 
and atom migration were investigated using the transition state search method. The results 
demonstrated that the CO on the uranium surface mainly manifests as chemical adsorption before 
dissociation of the CO molecule. The CO molecule exhibited a tendency to rotate and tilt upright 
adsorption. However, it is difficult for CO adsorption on the surface in one of the configurations 
with CO molecule in vertical direction but oxygen (О) 15 closer to the surface. Bader charge 
illustrates that the charge transfers from slab atoms to the 2л* antibonding orbital of CO molecule 
and particularly occurs in carbon (C) atoms. The time is less than 100 fs for the adsorptions that 
forms embryos with tilt upright in dynamics evolution. The density of states elucidates that the 
overlapping hybridization of C and О 2p orbitals is mainly formed via the d orbitals of uranium 
and molybdenum (Mo) atoms in the dissociation and re-adsorption of CO molecule. In conclusion, 
Mo-doping of the surface can decelerate the adsorption and dissociation of CO molecules. A 
Mo-doped surface, created through ion injection, enhanced the resistance to uranium-induced 
surface corrosion. 
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Introduction 


Uranium (U) holds the distinction of being the naturally occurring element with the highest 
atomic number (92) and relative atomic mass that has been discovered! It was found that 
uranium atoms release a large amount of energy through fission reactions. Uranium has been used 
in a wide range of industrial fields, such as nuclear power and aerospace""^!. There are three lattice 
structures of uraniumP!: alpha-uranium (a-U) with orthorhombic structure, beta-uranium (B-U) 
with body-centered tetragonal, and gamma-uranium (y-U) with body-centered cubic. The 
conversion temperatures from а to p, and В to y are 941 K and 1047 K, respectively. y-U lattice 
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structures are stabilized in uranium alloys by adding alloying elements with similar crystal 
structure and atomic size, including Molybdenum (Mo), zirconium (Zr), niobium (Nb), and 
titanium (Т1). U-Mo alloys are widely used in the nuclear industry because of their excellent 
corrosion resistance and mechanical properties. Therefore, it is considered a promising nuclear 
fuels/.10-11]. 

Surface corrosion is a typical problem owing to its special electronic structure and strong 
correlation with uranium 5f electrons in practical applications and the storage of 
uranium-molybdenum (U-Mo) alloys!?!, The electrochemical properties of uranium are as active 
as those of aluminum, both of which are susceptible to corrosion. Uranium can react with oxygen 
(O2), hydrogen (H2), water (H20), carbon dioxide (CO2), and carbon monoxide (CO) molecules at 
room temperature as per previous studies!?5, and hydrogen embrittlement can easily induce 
material powdering and burning due to the reaction between uranium and hydrogen!“ 
Additionally, uranium oxide is formed after the reaction between uranium and hydrogen. The 
oxide layer peels off after reaching a certain thickness, leading to further inward corrosion". 
Several researchers investigated the surface corrosion of uranium and uranium alloys. In the 
studies, reaction results were obtained with respect to small molecules interacting with material 
surfaces across various scales, from macro to micro, by utilizing techniques, such as Auger 


electron spectroscopy (AES)[18], and X-ray photoelectron spectroscopy (XPS) )!!?l, 


among 
others!?°l, Microscopic mechanisms were explored through surface reactions. For instance, 
Tellurium (Te) and О belong to the same group on the periodic table. A study by Wang et al. P! 
examined the effect of pure Nickel (N1) on the corrosion behavior of a GH3535 alloy in a Te vapor 
environment. Their findings suggested that the GH3535 alloy can resist Te corrosion, but at the 
expense of pure Ni. However, the study of surface corrosion using experimental methods has 
certain limitations due to localized experimental conditions and the unique properties of uranium. 

Recently, researchers have examined the theoretical mechanism of reactions between metals 
and small molecules. The interaction between hydrogen and 4d alloying atoms in Ni-based alloys 
has been studied using the density functional theory (DFT)P?. The larger size of the 4d atom 
induces stronger extrusion in the supercell, resulting in a greater repulsive force on the H atom. 
For uranium, researchers intend to study small-molecule adsorption and dissociation on uranium 
and U alloys using a theoretical method, which is considered the initial stage of surface corrosion. 
Research has been conducted on the performance and mechanisms of adsorption, dissociation, and 
diffusion, particularly with molecules such as O2, H2, and H20, on the surfaces of a-U, y-U, 
Uranium-titanium (U-Ti), and U-MoP??"l, The pronounced impact is attributed to the adsorption 
of carbon oxide on uranium and its alloys. 

The adsorption of carbon oxide on a U surface has not been as extensively studied as that of 
hydrogen and oxygen, despite attracting substantial theoretical interest. Prior studies primarily 
focused on static adsorption behavior, which involved determining different adsorption site 
configurations, studying structural relaxation, exploring adsorption energy at various sites, 
identifying optimal adsorption sites, and analyzing related parameters such as charge transfer, 
work function, and density of states. However, the dynamic behavior of small-molecule adsorption, 
molecule dissociation, and atom migration remain under-studied in extant studies, making it 
challenging to ascertain a comprehensive mechanism for the initial stage of surface corrosion. 
Consequently, it is crucial to consider the evolution of adsorption and dissociation in further 
studies. 


The evolution of the dynamics of CO adsorption and its dissociation on y-U (1 0 0) surface 
with different Mo doping was investigated using density functional theory, ab-initio molecular 
dynamics (AIMD), and climbing image-nudged elastic band method (CI-NEB)P?, This study is 
divided into four stages in this work: (1) Static calculations were performed on the adsorption of 
CO molecules in three reconstruction configurations for the Mo-doped U (1 0 0) surface, 
examining aspects, such as adsorption energy, configuration, bonding type, at all sites and optimal 
adsorption sites in different slabs. This foundational research is critical for subsequent analyses. (11) 
The study separately elucidates the crystal orbital Hamilton population (COHP), charge transfer, 
and total density of states (TDOS) of CO molecules at the optimal sites. (111) The evolution of the 
structure and two-dimensional differential charge density in the optimal adsorption models for CO 
adsorption is depicted using AIMD at 298.15 K. (iv) The minimum energy paths (MEPs) for CO 


molecule dissociation were explored. 


Numerical simulation details 


The Vienna ab-initio simulation package (VASP 5.4.4) was used for DFT, AIMD, and 
CI-NEB calculationsU???l, The projector-augmented wave (PAW) pseudopotential method is 
described in the electron-ion interaction section, suggesting an efficient plane-wave 
pseudopotential package in the electronic structure calculations of the surface model. The 
exchange and correlation energies of the electrons are determined by the Generalized Gradient 
Approximation (GGA) based on the Perdew (Burke) Ernzerhof exchange-correlation 
potential!?!**!_ This effect was neglected for spin polarization®*!. The plane-wave cut-off was set 
to 500 eV. 

The lattice parameter of y-U is optimized using a body-centered cubic (bcc) crystal cell, with 
an 11х11х11 Monkhorst-Pack k-point mesh employed for Brillouin zone integration. The lattice 
parameter is calculated to be 3.433 A, which aligns well with the theoretical value of 3.532 AM. 
A slab comprised of five U atomic layers is created to simulate the y-U (1 0 0) surface based on 
the computed lattice constant of y-U. The dimensions of the simulation box are 6.8668 A, 6.8668 
A, and 21.8668 A for length, width, and height, respectively. Three Mo doping configurations 
were constructed, namely: a single Mo atom replacing a U atom on the surface; a Mo atom 
substituting a U atom on the subsurface; and four Mo atoms replacing the entire top layer of U 
atoms on the surface. These reconstruction configurations were designated as Mo-U, Sec-Mo-U, 
and 4Mo-U for ease of subsequent analysis. 

The bottom three layers were fixed and top two layers were relaxed in the slab calculation for 
structural optimization апа AIMD simulations. K-point meshes of 5x5x1 and 7x7x1 were 
established for structural optimization and density of states, respectively. Additionally, the CO 
molecule was optimized and enclosed in a box with a side length of 10 A. The computed bond 
length was found to be 1.144 A, which closely aligns with the theoretical value of 1.129 A. A 
vacuum region along the z-axis was chosen to isolate the periodic images, maintaining a thickness 
of 15 A. Dipole correction was implemented to negate the dipole moment induced during CO 
adsorptionP?l, 


Adsorption energy (Fads) is obtained from the following equation: 
Es Esa Ea En (1) 


where Esia»+a denotes the total energy after CO adsorption and ЕА denotes the CO molecule energy. 


The negative adsorption energy demonstrates the stability of the adsorption system, and a larger 
absolute value of the adsorption energy indicates a more stable system. Conversely, a positive 
adsorption energy indicates that the system 15 unstable after adsorption. 

AIMD calculations were conducted using the Verlet algorithm over 1000 fs with a timestep 
of 1 fs, based on the canonical ensemble (NVT ensemble). The initial velocities of the CO 
molecules were assigned in accordance with the Maxwell-Boltzmann distribution at T — 298.15 K. 
A Langevin thermostatP9?7! was used to maintain the surrounding temperatures at 298.15 К, both 
at the start and end. The CI-NEB method was deployed to explore the diffusion of C and O atoms 
on the slab following CO adsorption. This method accurately determines the minimum energy 
paths and energy barrier, and also ensures the saddle point following convergence is the transition 
state. Energy convergence was established when the energy difference was less than 1x10 eV for 
adjacent iteration steps, and force convergence was defined as a residual force per atom less than 
2x10? eV/À. Particularly, CI-NEB calculations were repeated until the energy difference between 
adjacent iteration steps dropped below 1x107 eV. 


Results and discussion 


3.1 Configuration of adsorption and adsorption energy 

The adsorption of CO molecules on Mo-U, Sec-Mo-U, and 4Mo-U at highly symmetrical 
sites was investigated using the following calculations: The adsorption sites consist of top, hollow, 
and bridge sites. The adsorption sites on the Mo and U surface atoms must be considered under 
different chemical circumstances. The initial configuration comprised horizontal (Hor) and 
vertical (Ver) adsorption at each adsorption site. There are two adsorption categories for CO 
molecules, in which the O or C atoms are closer to the surface in the vertical direction. 

The top and side views of adsorption models of CO molecule on Mo-U, Sec-Mo-U and 
4Mo-U are shown in Fig.S1, Fig.S3 and Fig.S5, respectively, before structural relaxation. The top 
and side views of adsorption models of CO molecule on Mo-U, Sec-Mo-U and 4Mo-U are shown 
in Fig.S2, Fig.S4 and Fig.S6, respectively, after structural relaxation. The simplifications of the 
different adsorption configurations are described; for example, T-Hor, H-Hor, and B-Hor represent 
the top horizontal, hollow horizontal, and bridge horizontal, respectively. Тмо-Ног and Tu-Hor 
represent the top horizontal configurations of Mo and U atoms, respectively. Additionally, B- Ver 
and B- Ver-O represent the adsorption sites of the vertical bridge, where C and O atoms are closer 
to the surface, respectively. Furthermore, C and O atoms bonded to the slab-surface atoms. The 
adsorption of CO molecules occurs via chemical adsorption, but CO does not dissociate after 
adsorption, apart from the Hu-Hor and B-Hor configurations of Sec-Mo-U. This indicates that the 
CO molecule tends to rotate and tilt upright with stable adsorption and is approximately vertical to 
the surface in the final state. The C atom was close to the surface in the vertical direction. The 
phenomenon is also found in the CO molecule adsorption on a-U (0 0 1) surface??!. However, CO 
does not bond to slab atoms in the vertical adsorption configuration when the O atom is close to 
the surface. In some configurations, the CO molecule is far from the slab after structural 
relaxation. 

Fig.l shows a schematic of the CO single-molecule orbitalsP??!, The Highest Occupied 
Molecular Orbital (HOMO) and Lowest Unoccupied Molecular Orbital (LUMO) of free state CO 
molecule are the 5с orbital and 2л* orbital, respectively, as shown in Fig.l. The 5с orbital 
electrons were mainly concentrated near the C atomic side. The 2z* orbital is parallel to the C-O 


bond axis, which easily overlaps with the d orbitals of the U and Mo atoms on the slab. The CO 
molecule interacts with the substrate atoms mainly by bonding the C atom to the atoms, according 
to the matching requirement of the orbital orientation. Therefore, it can be inferred that the CO 
molecule is vertically adsorbed onto the surface in most configurations. 
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Fig.1. Schematic diagram of CO molecular orbitals 

Fig.2 shows the top and side views of the optimal adsorption configurations of CO molecules 
on the slab after structural relaxation. H-Hor (hollow horizontal), B-Hor (bridge horizontal), and 
H-Ver (hollow vertical) were the optimal adsorption configurations for Mo-U, Sec-Mo-U, and 
4Mo-U, respectively. The adsorption of CO molecules on slabs involves chemical adsorption, 
which is followed by bonding to the surface atoms. The adsorbed CO molecules tended to tilt 
upright on the surface in the optimal configurations of Mo-U and 4Mo-U. The CO molecule 
completely dissociates in the optimal configuration of Sec-Mo-U, and C atoms are bonded not 
only to surface atoms but also to Mo atoms on the subsurface. 

Tables S1—S3 show the adsorption energies and geometrical parameters of CO molecules оп 
Mo-U, Sec-Mo-U, and 4Mo-U, respectively. The chemical adsorption of the CO molecules 1s 
shown in Tables S1—S3. This indicates that dc.o increases when compared with the bond length of 
the CO-free state after adsorption, leading to the deterioration of the C—O bond strength. 

Table 1 lists the adsorption energy and geometrical parameters of the optimal configuration 
of CO molecules on the slab. The absolute value of the adsorption energy of Sec-Mo-U was larger 
than those of Mo-U and 4Mo-U, indicating that dissociated adsorption was more stable than 
undissociated adsorption. The absolute value of the adsorption energy of Mo-U was larger than 
that of 4Mo-U, indicating a decrease in Mo surface doping, leading to more stable CO adsorption. 
It is determined that the dpona between different bonding atoms in Table 1 is close to the bond 
length of U-C in UC (2.480 А)“, Mo-C in Мо›С (2.112 A)", and O-U in UO» (2.370 A), 
indicating that a stable chemical bond is formed between C or O atoms and U or Mo atoms on the 
slab. 


Fig.2. Top and side views of the optimal adsorption configurations of CO molecules оп the slab 
after structural relaxation: (a) Mo-U; (b) Sec-Mo-U; (c) 4Mo-U. 

Table 1 Adsorption energy and geometrical parameters in the optimal adsorption configurations of 
CO molecule adsorption on the slab: Faas (the adsorption energy), dc-o (the distance between C 
and О atom), Bond (the C and О atom bonding with the slab surface atom), свопа (bond length), 

and hc-s/ho-s (the distance between ће C/O atom and slab surface after adsorption) 


Configuration Exe V dc-o/ À Bond onal Å hc.s/À ho-s/À 
Mo-U/H-Hor -4.511 1.477 C-Mo 2.139 0.562 1.484 
C-U3 2.420 
C-Us 2.496 
O-U; 2.308 
O-U13 2.277 
Sec-Mo-U/B-Hor -6.146 3.905 C-Mo 2.264 0.383 0.595 
C-U3 2.406 
C-Us 2.341 
С-0 2.406 
C-Uio 2.341 
O-U; 2.213 
О-О 2.213 
4Mo-U/H-Ver -1.535 1.231 C-Moi 2.457 0.860 2.090 
C-Mo» 2.457 
C-Mo3 2.457 
C-Mo4 2.457 
C-U4 2.416 


Fig.3 depicts the COHP of the C and O atoms on the surface atoms for the optimal adsorption 
configurations of Mo-U, Sec-Mo-U, and 4Mo-U. The LOBSTER program"?! was used to 
calculate the COHP for the interatomic interactions. Positive and negative COHP results signify 
bonding and antibonding interactions, respectively. In Fig.3, the COHP of the C or O atoms is 
positive for the bonding atoms of the slab below the Fermi level, indicating bonding interactions 
between the CO molecule and slab atoms following CO adsorption. The strength of the interaction 
between two atoms is directly determined by the Integrated Crystal Orbital Hamilton Population 
(ICOHP)'5l, which is calculated using different energies below the Fermi level. An increase in 
bonding strength is associated with an increased ICOHP. Tables S3 and S4 display the ICOHP 
values for the COHP. As Table S4 shows, the ICOHP values for the CO molecule bonding to slab 
atoms follow a descending order of Sec-Mo-U, Mo-U, and 4Mo-U. This trend verifies the 
variation in adsorption stability of the optimal configuration in different slab models according to 
the ICOHP results. Additionally, Table S4 lists the ICOHP of C-O for the free-state and optimal 


adsorption configurations. The adsorption of CO molecules results in a reduction of the ICOHP of 
C-O. It was also observed that a higher interaction intensity between the C or O atoms and the slab 
atoms correlates with a larger degree of ICOHP degradation. 
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Fig.3. COHP of C and O atoms to the surface atoms which bond to CO molecule in the optimal 
adsorption configurations: (a) Mo-U; (b) Sec-Mo-U; (c) 4Mo-U 
Table S5-Table S7 shows the net charge distribution numbers for CO adsorption on Mo-U, 
Sec-Mo-U, and 4Mo-U for all adsorption configurations. The net charge is expressed as 
follows": 


datom M Bader Dvatence (2) 


where qaom, YBader, and valence denote the net atomic, Bader, and valence electrons, respectively. 
This indicated that the net charge of the CO molecule increased after adsorption, according to the 
charge-transfer results. Simultaneously, the slab lost electrons according to the principle of charge 
conservation, demonstrating that the electrons were mainly transferred from the slab atoms to the 
CO molecule. There are differences in the layers for the change in electrons after adsorption, 
causing the loss of the net charge number to mainly occur in the top two layers, especially in the 
first layer. It is believed that the C or O atoms mostly interact and transfer charge with the U and 
Mo atoms at the surface and subsurface. 

Table 2 lists the net charge distribution numbers for the optimal configurations of CO 
adsorption on the slab. The sum net charge number ranges from 0.9178 e to 2.7506 e for C and O 
atoms in the optimal models. The accumulation of charge is found to be closer to or larger than | e, 
indicating that the CO molecules form stable bonds with the slab atoms. The sum net charge 
number for the CO molecules exhibits a descending order with Sec-Mo-U, Mo-U, and 4Mo-U, 
which directly illustrates the interaction stability of the CO molecules on different slabs. 
Furthermore, the net charge gained by the C atoms is significantly larger than that by the O atoms, 
indicating that electrons are mainly transferred through the C atoms from the slab to the CO 
molecule. 

Fig. 4 displays the TDOS of the slabs both before and after CO adsorption in the most stable 
configuration. Upon CO adsorption on slabs, new peaks locally form within the ranges of -25 to 
-20 eV and -10 to -5 eV. Fig. 5 illustrates the TDOS of the CO molecule before and after its 
adsorption on the slab. It is observed that in the energy range of -5 to 5 eV, the HOMO and LUMO 
of the CO molecule exhibit significant interaction with slab atoms during CO molecule adsorption. 
The heights of the HOMO and LUMO orbitals undergo a substantial decrease, indicative of orbital 
delocalization. The Fermi level is elevated to a higher energy level, and the LUMO orbitals are 
filled with electrons due to electron transfer from the slab to the CO molecule. 

The integral value of the electron state in Fig.5 (c) is larger than those of (a) and (b) at energy 
levels ranging from -5 to 5 eV, leading to more nonlocalized electrons in the adsorption. The 


LUMO of CO molecule is 2л* orbital, which is comprised of ће p atomic orbital for С and О. It 
can be determined that the energy level position of 2л* orbital 15 closer to the C atomic orbital 

than that of О atomic in Fig.1. Therefore, the injecting electrons of 2л* antibonding orbital are 

mainly concentrated on the surrounding of C atoms, illustrating that the gaining of C atoms is 
more than that of O atoms in CO adsorption. 

Table 2 Net charge distribution number in the optimal configurations of CO adsorption on the slab: 
qc and qo denote net charge number of C апа О atoms, qtotaı represents the sum net charge number 
of C and О atoms, 415, q2nd, Q3:4, Q4, and qsa represent the total net charge number of first to fifth 

layers on the slab, respectively. 


Configuration qcle qole Quotal/e qis/e 4һа/е qxale а/е qsw/e 
Atom -1.0568 1.0568 0 - - - - - 
Mo-U - - - -0.0931 0.2681 -0.4432 0.8123 -0.5880 

Mo-U/H-Hor 0.6505 1.0540 1.7045 -1.7644 0.1914 -0.3812 0.8905 -0.6839 

Sec-Mo-U - - - -0.9969 1.3522 -0.5783 0.8912 -0.7120 
Sec-Mo-U/B-Hor 1.5348 1.2158 2.7506 -3.1005 0.5711 -0.5068 0.9115 -0.6704 
4Mo-U - - - 1.0060 -0.6764 -0.6701 0.8737 -0.5745 


4Mo-U/H-Ver -0.1328 1.0506 0.9178 0.2336 -0.8401 -0.6705 0.8444 -0.5256 
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Fig.4. TDOS of slabs before and after CO adsorption: before adsorption (a) Mo-U; (b) Sec-Mo-U; 
(c) 4Mo-U; after CO adsorption (a') Mo-U; (Б) Sec-Mo-U; and (с) 4Mo-U. 
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Fig.5. TDOS of CO molecule before and after adsorption on the slab: (a) before adsorption; (b) 
Mo-U; (c) Sec-Mo-U; (d) 4Mo-U 
3.2 Structural evolution of the optimal configurations 

Fig.6—Fig.8 present the structural evolution from top and side perspectives, as well as the 
two-dimensional differential charge density! in the optimal configurations of CO adsorption on 
Mo-U, Sec-Mo-U, and 4Mo-U, using AIMD simulation at 298.15 К. Concurrently, Figs. 57—89 
illustrate the corresponding three-dimensional differential charge density evolution at 298.15 K. 
CO progressively approaches the surface and uprightly rotates and tilts on Mo-U at 63 fs, resulting 
in a vertical adsorption stage of the CO molecule, where the C atom is closer to the surface. An 
approximate adsorption stage for 4Mo-U 15 observed at 83 fs, as shown in Fig.8. The CO molecule 
and slab atoms are respectively enveloped by red and blue regions, resulting in substantial electron 
accumulation and depletion. Charge accumulation is observed at C and O atoms enveloped by red 
regions. Conversely, charge depletion occurs at the slab-surface atoms adjacent to the C and O 
atoms, enveloped by blue regions. Additionally, at 0 fs in Figs. 6 (a)-8 (a), the interaction between 
the CO molecule and slab atoms is relatively weak. However, significant electron accumulation 
and depletion between the CO molecule and adjacent slab atoms can be observed as reaction time 
increases. After adsorption on Sec-Mo-U, the CO molecule is completely dissociated, as depicted 
in Fig.7, leading to the formation of bonds between the C and O atoms and the slab atoms. The 
spacing between C and O is 2.167 À at 48 fs, which suggests a break in the C-O bond, indicating 
a transition state between the molecular state and the dissociation of the CO molecule. The 
spacings between the C and O atoms further expand to 3.620 À and 4.239 À at 128 fs and 345 fs, 
respectively. The O and C atoms are eventually adsorbed at the hollow sites of the slab, in light of 
the rearrangement of the slab surface atoms. The differential charge densities of the C and O atom 
accumulation regions are higher than those of Mo-U and 4Mo-U in the final state, suggesting that 
a larger charge transfer enhances the stability of the C-U, O-U, and C—Mo bonds. 

Fig.9 presents the partial density of states (PDOS) at various evolution stages in the optimal 
configurations for CO adsorption at 298.15 K for (a) Mo-U, (b) Sec-Mo-U, and (c) 4Mo-U. The 
PDOS of the U and Mo atoms represents the PDOS of the surface and subsurface atoms bonded to 
C or O atoms in the slab, as depicted in Fig.9. The PDOS of U and Mo atoms extend beyond the 
Fermi level, indicating the strong metallic nature of U and Mo atoms. The PDOS hybridization 
occurs between C and O atoms, resulting in a distinct peak in the energy range from -12 to -4 eV 
at 0 fs. The PDOS of the C and O atom peaks diminish, and the peak broadens and overlaps with 
the U and Mo atoms. The peak values decrease following the broadening of the peaks for the 
PDOS of C and O atoms caused by the adsorption of CO molecules. The hybridization with U and 
Mo atoms highlights the bonding characteristics. Additionally, as shown in Fig.9 (a) and (c), no 
C-O bond breakage occurs after CO adsorption, followed by severe overlapping hybridization 
observable between the C and O atoms on Mo-U and 4Mo-U from the initial to the final state. 
However, at 48 fs, there is no evident hybridization of the PDOS between the C and O atoms, 
implying a weakening interaction between the C and О atoms. Eventually, after 122 fs, the 
hybridization completely vanishes between the C and O atoms, signifying a lack of interaction 
between C and O atoms. Independent peaks of the PDOS of the C and O atoms hybridize with the 
U and Mo atoms, respectively. 
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Fig.6 Evolution of structure with top and side views, and two-dimensional differential charge 
density in the optimal configuration of CO adsorption on Mo-U at 298.15 K. Red and blue color 
represent electrons accumulation and depletion, respectively. 
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Fig.7 Evolution of structure with top and side views, and two-dimensional differential charge 


density in the optimal configuration of CO adsorption on Sec-Mo-U at 298.15 К. Red and blue 


color represent electrons accumulation and depletion, respectively. 
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Fig.8 Evolution of structure with top and side views, and two-dimensional differential charge 
density in the optimal configuration of CO adsorption on 4Mo-U at 298.15 K. Red and blue color 
represent electrons accumulation and depletion, respectively. 
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Fig.9 PDOS of different evolution stages in the optimal configurations of CO adsorption at 298.15 
K: (a) Mo-U; (b) Sec-Mo-U; and (c) 4Mo-U 

3.3 Dissociation of CO molecule 

The minimum energy paths were calculated for CO molecule dissociation in the optimal 
configurations of Mo-U and 4Mo-U using the CI-NEB method, which provided a more profound 
understanding of the further dissociation mechanism of CO molecules on the U-Mo slab after 
adsorption. MEP values were calculated for the optimal configurations of Mo-U and 4Mo-U. 
Studies on the dissociation process of Sec-Mo-U are not discussed in this section because of the 
dissociation of CO molecule into C and O atoms on Sec-Mo-U in the optimal configuration. 
Transient calculations were based on the steepest descent method. Simultaneously, a specified 
number of images were linearly interpolated between the initial and final configurations. More 


importantly, spring forces were added to the constrained optimization of the images. 

Fig.10 and Fig.11 show the energy profiles for the dissociation of the CO molecule and the 
diffusion of atoms on Mo-U and 4Mo-U, respectively, in conjunction with the PDOS for the CO 
molecule and bonded atoms on Mo-U and 4Mo-U in various states. A potential barrier of 0.074 eV 
exists when the O atom overcomes Mo-U, leading to the rupture of the C-O bond. The potential 
barrier increases to 1.411 eV when the O atom surpasses the 4Mo-U. The O atoms migrate to 
adjacent sites, resulting in the bonding of the O atoms and slab surface atoms on Mo-U and 
4Mo-U in the final state. It is also noteworthy that the energy needed to break the C—O bond of the 
CO molecule adsorbed on Mo-U is less than that on 4Mo-U. This can be attributed to the weaker 
interaction of the C-O bond on Mo-U than on 4Mo-U after adsorption, as reflected by the ICOHP 
values for the C and O atoms on Mo-U and 4Mo-U in Table S3. 

The C 2p and О 2p orbitals hybridize with U 6p, U 6d, U 5f, Mo 4p, and Mo 4d orbitals in 
the region of -5 eV in initial state, as shown in Fig.10 (b) and Fig.11 (b). However, this peak was 
split, and its shape broadened after a decrease in the peak value. This suggested the formation of 
new chemical bonds in the transition state. There was a severe decrease in the PDOS overlapping 
hybridization of the C and O orbitals, implying that the interaction of C and O atoms became weak 
in the transient and final states. There was no obvious overlapping hybridization for the PDOS of 
C and O atoms in the final state, illustrating that the CO molecule was completely dissociated on 
Mo-U and 4Mo-U. Meanwhile, the PDOS of C and O atoms overlapping hybridize with U and Mo 
atoms separately, according to the overlapping of density of states in black or red dotted bordered 
rectangle in the ‘Final state’ in Fig.10 (b) and Fig.11 (b). Overlapping hybridization is mainly 
caused by the 2p orbitals of C and O atoms and d orbitals of U and Mo atoms during the 
dissociation and re-adsorption of CO molecules. 

Fig.S10 shows the energy of the system and maximum force of each atom versus the number 
of iteration steps in the transition-state calculation. This was terminated after 150 and 76 steps for 
the iterative calculations of Mo-U and 4Mo-U, respectively. The maximum force of each atom 
generally decreased with an increasing number of ionic steps. The energy of the systems gradually 
stabilized with an increase in the number of ionic steps for Mo-U and 4Mo-U. Both the force and 
energy reached the preset convergence standard, indicating that the iteration calculations were 
accurate for the CO molecule dissociation transition state on Mo-U and 4Mo-U. 
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Fig. 10. (a) Energy profile for CO molecule dissociation and atom diffusion on Mo-U and (b) 
PDOS of CO molecule and bonding atoms on Mo-U at different states 
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Fig.11. (a) Energy profile for CO molecule dissociation and atom diffusion on 4Mo-U (b) PDOS 
for CO molecule and bonding atoms on 4Mo-U at different states 


Conclusion 


The static and dynamic behaviors of CO adsorption on y-U (1 0 0) surface at different levels 
of Mo-doping were explored using DFT and AIMD. The adsorption of CO molecules on the slabs 
is predominantly chemical, except for the vertical orientation of the CO molecule. However, O 
atom is closer to the surface in the initial configuration. On these slabs, the adsorbed CO 
molecules do not tend to dissociate, but the CO molecules dissociate in the Hu-Hor and B-Hor 
configurations of Sec-Mo-U. Moreover, the optimal adsorption configurations were hollow 
horizontal, horizontal bridge, and hollow vertical, with adsorption energies of -4.511, -6.146, and 
-1.535 eV for Mo-U, Sec-Mo-U, and 4Mo-U, respectively. Bader charge analysis reveals that the 
charge transfer from surface atoms to the C atom exceeds that to the O atom, as the p orbital of the 
C atom is closer to the 2z* orbital than the p orbital of ће О atom. When the time was less than 
100 fs, the CO molecule on the slab demonstrated a nascent form of adsorption, leading to the 
vertical adsorption of CO molecules, where the C atom was closer to the surface. At 48 fs, in the 
B-Hor configuration of Sec-Mo-U, the spacing between C and О is 2.167 А, leading to the 
breakage of the C—O bond and dissociation of the CO molecule. 

There are two potential barriers of 0.074 eV and 1.411 eV encountered when the O atom 
surpasses Mo-U and 4Mo-U, respectively, leading to the breakage of the C-O bond. The O atoms 
migrate to the top of the Mo subsurface atoms and form a bridge between two adjacent Mo atoms 
on both Mo-U and 4Mo-U. Throughout the dissociation and re-adsorption of CO molecules, the 
2p orbitals of C and О atoms predominantly overlap with the d orbitals of U and Mo atoms. An 
increase in surface Mo doping not only elevates the difficulty of CO molecule adsorption but also 
raises the energy potential barrier of atom diffusion, thereby hindering the formation of an 
oxidation layer caused by CO adsorption on the uranium surface. This suggests that modulating 
Mo doping on the uranium surface can diminish the adsorption and dissociation of small 


molecules and enhance corrosion resistance during practical uranium storage and application. 
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Fig.S1. Top and side views of absorption models of CO molecule on Mo-U before structure 
relaxation: (a) Тмо-Ног; (b) Two- Ver; (c) Two-Ver-O; (d) Tu-Hor; (e) Tu-Ver; (f)Tu-Ver-O; (g) 
H-Hor; (h) H- Ver; (1) H-Ver-O; (j) B-Hor; (k) B-Ver; and (1) B-Ver-O. Blue, green, orange, and red 
balls represent the uranium, molybdenum, carbon, and oxygen atoms, respectively. 
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Fig.S2. Top and side views of absorption models of CO molecule on Mo-U after structural 
relaxation: (a) Тмо-Ног; (b) Тмо-Мет; (c) Two-Ver-O; (d) Tu-Hor; (e) Tu-Ver; (f)Tu-Ver-O; (g) 
H-Hor; (h) H-Ver; (i) H-Ver-O; (j) B-Hor; (К) B-Ver; and (1) B-Ver-O. Blue, green, orange, and red 
balls represent the uranium, molybdenum, carbon, and oxygen atoms, respectively. 

DE 


FEFFE: 


Fig.S3. Top and side views of absorption models of CO molecule on Sec-Mo-U before structural 
relaxation: (a) T-Hor; (b) T-Ver; (c) T-Ver-O; (d) Нмо-Ног; (e) Hmo-Ver; (f) Hmo-Ver-O; (g) Hu-Hor; 
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(b) Hu-Ver; (1) Hu-Ver-O; (j) B-Hor; (К) B-Ver; and (1) B- Ver-O. Blue, green, orange, and red balls 
represent the uranium, molybdenum, carbon, and oxygen atoms, respectively. 


Fig.S4. Top and side views of absorption models of CO molecule on Sec-Mo-U after structural 
relaxation: (a) T-Hor; (b) T-Ver; (c) T-Ver-O; (d) Нмо-Ног; (e) Hmo- Ver; (f) Hwo-Ver-O; (g) Hu-Hor; 
(h) Hu-Ver; (1) Hu-Ver-O; (j) B-Hor; (К) B-Ver; and (1) B-Ver-O. Blue, green, orange, and red balls 

represent the uranium, molybdenum, carbon, and oxygen atoms, respectively. 
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Fig.S5. Top and side views of absorption models of CO molecule on 4Mo-U before structural 
relaxation: (a) T-Hor; (b) T-Ver; (c) T-Ver-O; (d) H-Hor; (e) H-Ver; (f) H-Ver-O; (g) B-Hor; (h) 
B- Ver; and (h) B-Ver-O. Blue, green, orange and red balls represent the uranium, molybdenum, 

carbon, and oxygen atoms, respectively. 
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Fig.S6. Top and side views of absorption models of CO molecule on 4Mo-U after structural 
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relaxation: (a) T-Hor; (b) T-Ver; (c) T-Ver-O; (d) H-Hor; (e) H-Ver; (f) H-Ver-O; (g) B-Hor; (h) 
B- Ver; and (h) B-Ver-O. Blue, green, orange and red balls represent the uranium, molybdenum, 
carbon, and oxygen atoms, respectively. 


iu Ф 02° а N VS or 


X5 IE 


Fig.S7. Top and side views of three-dimensional differential charge density evolution in the most 
stable models of CO adsorption on Mo-U at 298.15 K: (a) 0 fs; (b) 63 fs; (c) 186 fs; and (d) 305 fs. 
Red and yellow indicate an increase and a decrease, respectively, in charge density (Isosurface 
level: 0.004 e/A?). Blue, green, orange, and red balls represent uranium, molybdenum, carbon, and 
oxygen atoms, respectively. 
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Fig.S8. Top and side views of three-dimensional differential charge density evolution in the most 
stable models of CO adsorption on Sec-Mo-U at 298.15 K: (a) 0 fs; (b) 48 fs; (c) 122 fs; and (d) 
358 fs. Red and yellow indicate an increase and a decrease, respectively, in charge density 
(Isosurface level: 0.004 e/À?). Blue, green, orange, and red balls represent the uranium, 


molybdenum, carbon, and oxygen atoms, respectively. 
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Fig.S9. Тор and side views of three-dimensional differential charge density evolution in the most 

stable models of CO adsorption on 4Mo-U at 298.15 K: (a) 0 fs; (b) 83 fs; (c) 331 fs; and (d) 680 

fs. Red denotes an increase in charge density (Isosurface level: 0.002 e/À?), and yellow indicates a 
decrease in charge density (Isosurface level: 0.004 e/À?). Blue, green, orange, and red balls 


represent uranium, molybdenum, carbon, and oxygen atoms, respectively. 
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Fig.S10. Energy of the system, and maximum force of each atom versus the number of iteration 
steps in transition state calculation: (a) Mo-U and (b) 4Mo-U 
Table S1 Adsorption energy and geometrical parameters of CO molecule adsorption on Mo-U: 
Eads ( adsorption energy), dc-o ( distance between C and О atoms), Bond ( C and O atoms bonding 
with the slab surface atom), dona (bond length), and Ac-s/ho-s ( distance between the C/O atom and 
the slab surface after adsorption) 


Configuration Ease V dc-olÅ Bond onal Å hc-s/Å hos/A 
Тмо-Уег -1.840 1.198 С-Мо 2.323 1.700 2.703 
C-Ui5 2.396 
Two- Ver -1.567 1.175 C-Mo 2.040 1.961 3.136 
Two-Ver-O 0.051 1.156 - - 3.456 2.299 
Tu-Hor -2.054 1.286 C-U3 2.364 1.709 2.311 
C-Us 2.194 
O-U; 2.261 
Tu- Ver -0.781 1.171 C-Us 2.387 1.864 3.035 
Tu-Ver-O -0.010 1.144 - - 4.416 3.270 
H-Hor -4.511 1.477 C-Mo 2.139 0.562 1.484 


H-Ver 


H-Ver-O 


B-Hor 


B-Ver 


B-Ver-O 


-3.969 


-0.966 
-4.444 


-2.271 


-2.431 


1.396 


1.145 
1.466 


1.244 


1.244 


O-U; 
O-Us 
C-Mo 
C-Us 
O-Mo 
C-Mo 
O-Us 


2.496 
2.308 
2.277 
2.126 
2.421 
2.409 


2.145 
2.425 
2.467 
2.265 
2.302 
2.175 
2.305 
2.339 
1.961 
2.445 


0.620 


4.402 
0.568 


1.417 


1.611 


1.768 


3.630 
1.495 


2.191 


2.298 


Table S2 Adsorption energy and geometrical parameters of CO adsorption on Sec-Mo-U 


Configuration Eaas/e V dc-o/ À Bond dBosa/ Å hic-s/À ho-s/Å 
T-Hor -4.120 1.509 C-Mo 2.325 0.573 1.410 
C-Us 2.403 
C-Un 2.274 
O-Us 2.247 
О-0 2.243 
T-Ver -2.888 1.375 C-Us 2.480 0.806 1.161 
C-Un 2.480 
C-Uio 2.226 
O-U3 2.296 
T-Ver-O -0.794 1.150 - - 4.087 2.953 
Hmo-Hor -4.309 1.477 C-Mo 2.321 0.543 1.401 
C-Us 2.395 
C-Us 2.500 
C-Un 2.395 
C-Uio 2.500 
O-Us 2.245 
O-Ui9 2.245 
Нмо-Мег -1.395 1.239 - - 1.342 2.581 
Нмо-Мег-О 0.031 1.146 - - 4.436 3.290 
Hu-Hor -5.716 4.299 C-Us 2.476 0.270 1.235 
C-Us 2.252 
C-Uio 2.252 
O-U3 2.111 
О-О» 2.111 
Hu-Ver -1.109 1.242 - - 1.294 2.536 
Hu-Ver-O -0.001 1.144 - - 4.893 3.749 
B-Hor -6.146 3.905 C-Mo 2.264 0.383 0.595 
C-Us 2.406 
C-Us 2.341 
C-Un 2.406 
C-Uio 2.341 
O-U3 2.213 
О-О» 2.213 
B-Ver -1.759 1.198 C-Us 2.473 1.792 2.984 
C-Uio 2.473 


B-Ver-O -0.546 1.144 - - 4.817 3.694 
Table S3 Adsorption energy and geometrical parameters of CO adsorption on 4Mo-U 


Configuration Exas/eV ас-о/А Вопа dona А hcs/À ho-s/Å 
T-Hor -1.466 1.164 C-Mo4 2.070 2.135 3.299 
T-Ver -1.464 1.164 C-Mo4 2.070 2.137 3.301 

T-Ver O -0.102 1.153 - - 3.580 2.427 
Н-Ног -1.532 1.231 C-Moi 2.483 0.864 2.093 
C-Mo? 2.437 
C-Mo3 2.483 
C-Mo4 2.437 
H-Ver -1.535 1.231 C-Moi 2.457 0.860 2.090 
C-Mo? 2.457 
C-Mo3 2.457 
C-Mo4 2.457 
C-U4 2416 
H-Ver-O 0.019 1.144 - - 4.410 3.266 
B-Hor -1.333 1.210 C-Mo? 2.199 1.189 2.399 
C-Mo4 2.199 
B-Ver -1.332 1.210 C-Mo? 2.198 1.181 2.391 
C-Mo4 2.198 
B-Ver-O -0.003 1.144 - - 4.537 3.393 


Table S4 ICOHP for C and O with surface atoms of the optimal adsorption configurations in 
Mo-U, Sec-Mo-U, and 4Mo-U 


Configuration Bond (C atom) ICOHP/eV Bond (O atom) ICOHP/eV 
Free state / / O-C 20.3919 
Mo-U C-U3 2.1566 O-C 9.3406 
С-05 3.0340 O-U3 2.1532 
C-Mo 3.3764 O-U13 5.7301 
Sec-Mo-U C-U3 5.4662 O-C 0.0007 
C-U8 4.6889 O-U3 7.2016 
C-U12 5.6177 O-U12 7.4664 
C-U19 5.1382 
C-Mo 2.4053 
4Mo-U C-Mol 1.2833 O-C 17.2081 
C-Mo2 1.4468 
C-Mo3 1.4468 
C-Mo4 1.6366 
C-U4 4.4582 


Table S5 Net charge distribution number of CO adsorption on Mo-U. Specifically, qc and qo are 
the net charge numbers of C and О atoms, respectively, qtotai is the sum of net charge numbers of C 
and О atoms, qist, q2nd, 3га, фа, and qsm denote the total net charge number of first to fifth layers 
on the slab, respectively. 


Configuration qcle qole total/e qis/e q2nd/€ дэге qatn/e qsa/e 
Atom -1.0568 1.0568 0 - - - - - 
free surface - - - -0.0931 0.2681 -0.4432 0.8123 -0.5880 


Тмо-Ног -0.4149 1.0723 0.6574 -0.8330 0.2740 -0.3173 0.7311 -0.5541 


Two- Ver -0.6307 1.1061 0.4754 -0.5503 0.2081 -0.4121 0.8184 -0.5816 
Тмо-Мег-О -0.9541 1.1039 0.1498 -0.1836 0.1890 -0.4310 0.8182 -0.5877 
То-Ног 0.0121 1.1053 1.1174 -1.2675 0.1950 -0.3063 0.8578 -0.6374 
Tu-Ver -0.6237 1.0487 0.4250 -0.2625 0.0870 -0.5326 0.8078 -0.5666 
Tu-Ver-O -1.1183 1.1519 0.0336 -0.1334 0.2750 -0.4473 0.8109 -0.5883 
H-Hor 0.6505 1.0540 1.7045 -1.7644 0.1914 -0.3812 0.8905 -0.6839 
H-Ver 0.4977 1.0629 1.5606 -1.5717 0.1378 -0.3977 0.8756 -0.6467 
H-Ver-O -1.1016 1.1237 0.0221 -0.3340 0.4746 -0.4110 0.8391 -0.6397 
B-Hor 0.6523 1.0529 1.7052 -1.7517 0.1679 -0.3803 0.8864 -0.6699 
B-Ver -0.1173 1.0683 0.9510 -1.0545 0.1498 -0.2799 0.7619 -0.5695 
B-Ver-O -0.2867 1.1373 0.8506 -1.1311 0.3151 -0.2546 -0.5842 -1.9324 
Table S6. Net charge distribution of CO adsorption on Sec-Mo-U 
Configuration qcle qole total/e qis/e 4һа/е qxa/e qane qsa/e 
Atom -1.0568 1.0568 0 - - - - - 
free surface - - - -0.9969 1.3522 -0.5783 0.8912 -0.7120 
T-Hor 0.8038 1.0272 1.8310 -2.3498 0.7630 -0.4986 0.8351 -0.6234 
T-Ver 0.4456 1.1390 1.5846 -2.0926 0.7752 -0.4857 0.8323 -0.6560 
T-Ver-O -1.0578 1.1376 0.0798 -0.9726 1.1777 -0.5241 0.8208 -0.6292 
Hwo-Hor 0.7371 1.0355 1.7726 -2.3345 0.8240 -0.5410 0.9187 -0.6823 
Hwo- Ver 0.0290 0.9949 1.0239 -1.6528 0.8909 -0.4859 0.8119 -0.6300 
Hwo- Ver-O -1.0770 1.1193 0.0423 -1.0348 1.3368 -0.5733 0.8923 -0.7121 
Hu-Hor 1.5379 1.1452 2.6831 -3.1190 0.7488 -0.5584 0.8857 -0.6830 
Hv- Ver 0.0126 1.0047 1.0173 -1.8713 1.2154 -0.5569 0.9043 -0.7504 
Hu-Ver-O -1.1044 1.1179 0.0135 -1.0151 1.3491 -0.5799 0.8939 -0.7131 
B-Hor 1.5348 1.2158 2.7506 -3.1005 0.5711 -0.5068 0.9115 -0.6704 
B-Ver -0.3564 1.0245 0.6681 -1.4956 1.1805 -0.6399 0.9229 -0.6779 
B-Ver-O -1.0738 1.0919 0.0181 -0.9294 1.2670 -0.6527 0.9038 -0.6584 
Table S7 Net charge distribution for CO adsorption on 4Mo-U. 

Configuration qcle qole Qtota/ € qis/e Q2»d/€ qale qa/e qsa/e 
Atom -1.0568 1.0568 0 - - - - - 
free surface - - - 1.0060 -0.6764 -0.6701 0.8737 -0.5745 

T-Hor -0.7455 1.0654 0.3199 0.6723 -0.6999 -0.6414 0.8867 -0.5798 
T-Ver -0.7635 1.0888 0.3253 0.6721 -0.7023 -0.6464 0.8882 -0.5794 
T-Ver-O -0.9979 1.0810 0.0831 0.9620 -0.7413 -0.6756 0.9057 -0.5806 
H-Hor -0.0710 0.9892 0.9182 0.2338 -0.8338 -0.6792 0.8479 -0.5272 
H-Ver -0.1328 1.0506 0.9178 0.2336 -0.8401 -0.6705 0.8444 -0.5256 
H-Ver-O -1.0818 1.0922 0.0104 0.9872 -0.6638 -0.7215 0.9123 -0.5749 
B-Hor -0.2931 0.9341 0.6410 0.5519 -0.8874 -0.6935 0.9475 -0.5998 
B-Ver -0.2966 0.9402 0.6436 0.5868 -0.9285 -0.6727 0.8984 -0.5679 
B-Ver-O -1.1184 1.1250 0.0066 1.0116 -0.6835 -0.7037 0.8942 -0.5763 


